PROCEEDINGS 

0F SCIENCE 



DESY 09-196 
Edinburgh 2009/20 
LTH 853 
TUM-T39-11-11 



Lattice Investigations of Nucleon Structure at Light 
Quark Masses 

ON 

o 
° 

<N M. Gockeler", Ph. Hagler /; , R. Horsley c , Y. Nakamura", D. Pleiter^, P. E. L. Rakow e , 

O A. Schafer G. Schierholz a -f, H. Stiiben^, J. M. Zanottr 

a Institut fiir Theoretische Physik, Universitat Regensburg, 93040 Regensburg, Germany 
h Institut fiir Theoretische Physik T39, TU Miinchen, 85747 Garching, Germany 
c School of Physics and Astronomy, University of Edinburgh, Edinburgh EH9 3JZ, UK 
i— i d John von Neumann Institute NIC / DESY Zeuthen, 15738 Zeuthen, Germany 

e Theoretical Physics Division, Department of Mathematical Sciences, University of Liverpool, 
T 1 Liverpool L69 3BX, UK 

Ph / Deutsches Elektronen-Synchrotron DESY, 22603 Hamburg, Germany 

pH s Konrad-Zuse-Zentrum fiir Informationstechnik Berlin, 14195 Berlin, Germany 

l— - 1 E-mail: jzanotti@ph.ed.ac.uk 

QCDSF/UKQCD Collaboration 

NO Lattice simulations of hadronic structure are now reaching a level where they are able to not only 

f-*>) complement, but also provide guidance to current and forthcoming experimental programmes at, 

^vq e.g. Jefferson Lab, COMPASS/CERN and FAIR/GSI. By considering new simulations at low 

quark masses and on large volumes, we review the recent progress that has been made in this 
exciting area by the QCDSF/UKQCD collaboration. In particular, results obtained close to the 
physical point for several quantities, including electromagnetic form factors and moments of ordi- 
JTj nary parton distribution functions, show some indication of approaching their phenomenological 

s_! values. 



The XXVII International Symposium on Lattice Field Theory - LAT2009 
July 26-31 2009 

Peking University, Beijing, China 
* Speaker. 



© Copyright owned by the author(s) under the terms of the Creative Commons Attribution-NonCommercial-ShareAlike Licence. 



http: //pos. sissa. it/ 



Nucleon Structure at Light Quark Masses 



J. M. Zanotti 



U 



1.6 



O 1.4 



< 



1.2 



1.0 



0.0 



▼ 



0.2 



0.4 

ml [GeV 2 ] 



0.6 



▼ 5.25 
A 5.29 
♦ 5.4 



Oi 



Figure 1: Nucleon mass as a function of nv\ for three different values of j3. The scale is set using r§ 
0.467 fm and the physical mass is indicated by the cross. 



1. Introduction 



Much of our knowledge about hadronic structure in terms of quark and gluon degrees of free- 
dom has been obtained from high energy scattering experiments. However, there are still many 
unresolved issues in hadronic physics that need to be addressed, from both an experimental and 
theoretical perspective. This is one of the main motivations of the 12 GeV Jefferson Lab upgrade 
which aims to [1]: search for exotic mesons; study the role of hidden flavours in the nucleon; map 
out the spin and flavour dependence of the valence partem distribution functions; explore nuclear 
medium effecs; and measure the generalised parton distribution functions of the nucleon. It is im- 
perative that these and other exciting experimental efforts, such as those at COMPASS/CERN and 
FAIR/GSI, are matched by modern lattice simulations. See recent lattice reviews [2, 3, 4, 5] for 
latest developments. 

In this talk we present the latest results from the QCDSF/UKQCD Collaboration for elec- 
tromagnetic (EM) form factors, (x) q , (x)Aq, gA and gj. Simulations are performed with the 
Wilson gauge action and two flavours of ^?(a)-improved Wilson fermions. Four values of /3 = 
5.20, 5.25 ,5.29 ,5.40, corresponding to lattice spacings in the range 0. 1 < a < 0.07 fm, allow for 
the approach to the continuum limit to be assessed, while a range of lattice volumes (1.1 < L < 
3.2 fm) enable us to search for finite size effects in our simulations. A key feature of the new results 
presented here is the control we are now able to achieve over the approach to the chiral limit. This 
is made possible by several new simulations with pion masses reaching as low as 170 MeV, which 
are summarised in Fig. 1. 
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Figure 2: Results for the isovector Pauli radius, r%. 



Figure 3: Isovector anomalous magnetic moment, K v 
in nuclear magnetons. 



2. Electromagnetic Form Factors 

The study of the electromagnetic properties of hadrons provides important insights into the 
non-perturbative structure of QCD. The EM form factors reveal important information on the in- 
ternal structure of hadrons including their size, charge distribution and magnetisation. 

A lattice calculation of the g 2 -dependence of hadronic electromagnetic form factors can not 
only allow for a comparison with experiment, but also help in the understanding of the asymptotic 
behaviour of these form factors, which is predicted from perturbative QCD. Such a lattice calcula- 
tion would also allow for the extraction of other phenomenologically interesting quantities such as 
charge radii and magnetic moments. For a recent review see [6]. 

On the lattice, we determine the form factors F\{q 2 ) and Fj{q 2 ) by calculating the following 
matrix element of the electromagnetic current 

{p',s'\nq)\p,s) =u{p\s') 

where u(p, s) is a Dirac spinor with momentum, p, and spin polarisation, s, q = p' — p is the 
momentum transfer, is the nucleon mass and is the electromagnetic current. The Dirac (F\ ) 
and Pauli [F%) form factors of the proton are obtained by using m' = \uy^u — hdy^d, while for 
isovector form factors f„ = uy^u — dy^d. It is common to rewrite the form factors F\ and F2 in 
terms of the electric and magnetic Sachs form factors, G e = F\ +q 2 / (2M^) 2 F2 and G m = F\ +F2. 

If one is using a conserved current, then (e.g. for the proton) F^ (0) = (0) = 1 gives the 
electric charge, while G$ (0) = pL ^ = 1 + fc^ gives the magnetic moment, where F^^ (0) = 
is the anomalous magnetic moment. From Eq. (2. 1) with see that F2 always appears with a factor 
of q, so it is not possible to extract a value for F2 at q 2 = directly from our lattice simulations. 
Hence we are required to extrapolate the results we obtain at finite q to q = 0. Form factor radii, 
r i = "J ( r i)' we defined from the slope of the form factor at q 2 = 0. 

In Fig. 2 we see results for the isovector Pauli radius plotted as a function of m 2 . Here we 
observe a common trend of lattice determinations of this quantity, namely that the results are a 
factor of ~ 3 smaller than experiment and show little variation as a function of m| for pion masses 



y»F l (q 2 )+io» v ^ N F 2 (q 2 ) u(p, s) , 



(2.1) 
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larger than m n > 300 MeV. While the new results below m n < 300 MeV are showing signs of 
upward curvature, as predicted from chiral perturbation theory [7, 8, 9], it is not yet clear whether 
this is enough to reach the experimental value. Of course, at such light quark masses, finite volume 
effects are expected to suppress the nucleon charge radii [9], and the results from three different 
volumes at m\ « 0.09 GeV 2 indicate that this is indeed the case, although the effect appears to be 
rather small at this pion mass. 

Results for the isovector anomalous magnetic moment, ?c v , are shown in Fig. 3 and the results 
are seen to be slowly increasing towards the experimental value as the pion mass is decreased. 
Closer to the chiral limit, the results are stalling to exhibit stronger chiral behaviour, once again as 
predicted by %PT [7, 8, 9]. Here finite volume effects are also expected to suppress the magnetic 
moment [8], so the fact that our results on a finite volume are not increasing fast enough to agree 
perfectly with experiment is not totally unexpected. 

3. Moments of Structure Functions 

In this section we discuss our latest findings for the nucleon's axial charge, gA, tensor charge, 
gr, and spin-independent and spin-dependent quark momentum fractions, (x) q and {x)^. All 
operators in this section have been renormalised nonperturbatively using the Rome-Southampton 
method [10]. 

3.1 Axial Charge, gA 

The axial coupling constant of the nucleon is important as it governs neutron j3 -decay and 
also provides a quantitative measure of spontaneous chiral symmetry breaking. It is also related to 
the first moment of the helicity dependent quark distribution functions, gA = Au — Ad. It has been 
studied theoretically as well as experimentally for many years and its value, gA = 1.2695(29), is 
known to very high accuracy. Hence it is an important quantity to study on the lattice, and since it 
is relatively clean to calculate (zero momentum, isovector), it serves as useful yardstick for lattice 
simulations of nucleon structure. 

The axial charge is defined as the value of the isovector axial form factor at zero momentum 
transfer and is determined by the forward matrix element 

(p,s\A u - d \p,s) =2g A s fl , (3.1) 

where = qj^Jsq, p is the nucleon momentum, and is a spin vector with s 2 = — M^. 

g A has been studied in-depth for many years by the QCDSF/UKQCD [11], LHP [12] and 
RBC/UKQCD [13] collaborations and has been shown to suffer from large finite size effects. Such 
effects are also seen in our new results at lighter pion masses as seen in Fig. 4 where we show 
those results for where we have more than one volume available at a fixed pion mass. Here we 
normalise the finite volume results for the axial charge, gA{L), with the result obtained from the 
largest available volume for that choice of (j3, k), gA{°°)- 

Another way to visualise our results is to consider the ratio gA/ fit, as shown in Fig. 5. A 
good reason for considering this ratio is that the common renormalisation constant, Za, cancels in 
the ratio. It is also a natural ratio to consider in the context of the Goldberger-Treiman relation, 
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Figure 4: Finite volume results for the nucleon ax- 
ial charge, £a(L), normalised with the largest avail- 
able volume, gA(°°), plotted as a function of m n L 
for J3 =5.29. 
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Figure 5: Results for the ratio of the nucleon ax- 
ial charge and the light pseudoscalar decay constant, 

gA/fn 



gA = fngnNN /Mn, which is valid in the chiral limit and approximately valid at larger quark masses. 
As can be seen in Fig. 5, the results for gA/fn show a strong dependence on m\ and increase 
towards the experimental value as they approach the physical point. This may be an indication 
that the leading finite size effects cancel in this ratio [11, 14]. It would be nice to see if such 
cancellations occur in #PT expressions of these quantities. This is currently under investigation. 

3.2 Tensor Charge, gr 

Although the tensor charge has not received as much attention as gA, there have been recently 
some new attempts to calculate the quantity on the lattice [15, 16, 17]. Since gr is not known 
experimentally, it provides an opportunity for the lattice to make a prediction, although given the 
difficulties with gA as discussed above, care must be taken, not only with the chiral extrapolation, 
but also in the assessment of finite size effects. 

In Fig. 6 we show the latest QCDSF/UKQCD results, including the new simulations at low 
pion masses. Here we observe a similar trend to that observed in the axial charge, gA, above, 
namely that increasing the volume at a fixed pion mass increases the result, although the effects are 
not as strong. 

Considering the behaviour of the results at fixed m n L, we observe that the results are reason- 
ably flat as a function of m 2 n . This is to be compared with the above result for gA where we observed 
a definite increase in the results. 

Given that we still don't have the chiral and finite volume systematics under control, it is not 
yet possible to make a strong prediction for the isovector tensor charge of the nucleon, however our 
current results indicate that it is likely to be in the range 0.9 < gr < 1.1. 

3.3 Nucleon Momentum Fraction, (x) 

Lattice studies of (x) q are notorious in that all lattice results to date at heavy quark masses 
exhibit an almost constant behaviour in quark mass towards the chiral limit and are almost a factor 
of 1.5 larger than phenomenologically accepted results, e.g. {x)^ R j T = 0.157(9) [18]. Despite 
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Figure 6: Results for the isovector tensor charge, #7-, Figure 7: Results for the isovector quark momentum 
in the MS scheme at /I = 2 GeV. fraction, (x) H _d, in the MS scheme at /I = 2 GeV. 



predictions that the results should drop dramatically in the light quark mass region [19], we are still 
yet to see conclusive evidence for this from a lattice calculation. 

To date, only connected contributions have been simulated to high precision, and although 
progress is being made towards the evaluation of disconnected contributions [20], here we only 
quote results for isovector quantities where disconnected contributions cancel. 

Figure 7 shows the latest QCDSF/UKQCD results including new results below m % < 300 MeV. 
While the results above nv\ > 0.1 GeV 2 show the same constant behaviour seen in many lattice 
simulations, the new results in the region below m 2 < 0.1 GeV 2 are beginning to show some inter- 
esting behaviour, albeit with larger statistical errors. What is particularly interesting is the results 
at a fixed pion mass (ps 270 MeV), but with three different lattice volumes (L = 24, 32, 40). Here 
we observe that the result from the smallest volume continues the trend of the results at heavier 
pion masses, while the results from the larger volumes sit slightly lower. This is not totally un- 
expected, as the dramatic downward curvature has been shown to be only achievable on larger 
volumes [9, 21, 22]. For this reason, we expect that even after a large increase in statistics, the 
current point at m % 170 MeV with L = 40 could well still lie somewhat higher than phenomeno- 
logical determinations. A further simulation at this pion mass but with spatial volume length of 
L = 64 is currently underway to further test these ideas. 

Similar behaviour is seen in the spin-dependent momentum fraction, (jc)a^. 



4. Conclusion & Outlook 

We have presented some of the latest results on nucleon structure from the QCDSF/UKQCD 
collaboration, including EM form factors, (x) q , (x)Aq, gA and gj. These calculations are now be- 
coming available at pion masses as low as m K rj 170 MeV, so direct comparison with experimental 
determinations will soon be possible. However, as we have seen in, e.g. gA, finite size effects are 
starting to become a serious issue. As a result, we are now planning a new simulation on a volume 
of (5fm) 3 , in order to minimise these effects, although corrections from ^PT will still probably 
need to be taken into account. 
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